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ABSTmCT 

The removal of phosphorus fron soils was shewn to be the 
result of three mechanisms: ion exdiange, physical adsorption and 
precipitation. The degree of sorption by the combination of these 
mechanisms was evaluated by an equilibriijii-adsorption procedure. 
TTie equilibrium distribution between the soil and liquid phases 
were correlated using the Lanptiuir Isotherm Equation. The Lan^ouir 
Parameters were incorporated in the expression 



V /v = 1 + 

0) c 



o 

which gave an estimate of the rate of phosphorus movement in the 
soil. The life-time of a soil -groundwater system was then determned 
from the groundwater flow rate and the calculated water-to-contaminant 
velocity ratio. 
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1.0 INTRODUCTION 



The treatment of water by adsorption techniques is an 

ancient process as indicated by records as far back as the time 

(1) . . 

of Aristotle. At this time, sand filters were used for the 

purification of sea water and impure drinking waters. According 
to Bacon in Sylva SylvaAum : 



". ....... . ,tAAjat hath been made, ojf ^aJtt MiteA. paJi6Q-d 

th/ioagh aahth thn^oagh ten ue64e^, one Mttkln anotheA, and yet 
tt hath not lo6t Xt4 6aZttne^6 a6 to becorm potabte., but whe.n 
d/iayned thxaugh tiA)anty fe^4e£4 hath bo^aamt {^KdMh," 



The ability of clays and soils to adsorb coraponents of manure 

liquors was studied extensively by Sir Hunphrey I!a\y, Lantiuschini , 

(2) 
Huxtable and others m the early part of the nineteenth century. 

Credit for the recognition of ion exchange phenomena is generally 

(3) 

attributed to two English agricultural diemists, TTicmpson and 

Way. 

Today, purification of septic tank effluent from private 
waste sources is accomplished chiefly by sand or soil percolation 
for sites beyond reach of municipal sewage facilities. In this 
application, the soil acts both as a filter and an ion exchange- 
adsorption medium. The efficiency of a septic tank tile field 
systen is heavily dependent on the capability of the soil to remove 
contaminants which may be detrimental to neighbouring sources of 
natural water. Thus, the movement of these contaminants throu^ 
the soil poses a nHJor difficulty in the design and planning of 
septic tank systens. 



2.0 ION EXCHANGE MD ADSORPTION 



2.1 Structui?e of Ion Exchangers 

The natural materials responsible for the ion=exchange 

phenomena have been identified as clays, glauconites, zeolites and 

(14) 
huno-C acids. According to Helffench, the most natural ion-exchange 

minerals are crystalline aluminosilicates with cation exchange pro- 
perties. Zeolites are typical representatives of this group of 
materials which include the minerals analcite, Na(Si^A10„)-H„05 

chabazite, (Na, Ca) , CSi^A10c)o'6HoO, harmotome, (K, B ) (SicA1^0T,.)-H„0 

/ b / / — a o Z 14- i 

heiilandite, Ca(Si_A10o)-5H^0, and natrolite, Na^(Si.A1^0, ^)*2H„0.: 
382 232 lU 2 

All these minerals have a relatively open thiree dimensional framework 
structure with diannels and interconnecting cavities in the alumino- 
silicate lattice. The zeolite lattice consists of SiO,, and AID, 
tetrahedra which have their o>c/gen atoms in conmon. Since aluminum 
is trivalent, the lattice carries a negative electric charge (one 
elerentary charge per aluminum atom) . This charge is balanced by 
alkali or alkaline earth cations whidi do not occupy fixed positions 
but are free to move in the channels of the lattice framework. These 
ions act as counter ions and can be replaced by other cations. 

Certain alundnosilicates can also act as anion exchangers, 

- 2- 3- 
The exchange of OH for CI , S0^ , and PO^ has been observed m 

mon'taiDrillonite J kaolinite, and in feldspars belonging to the 

sodalite and cancrinite groups. 

In general, clays originate as components in residual 

bodies through weathering of a geological feat\ire or frxm hydrx^thermil 

synthesis at higl'i tenperature cand pressure. A detailed discussion of 

. (8) 
"tiiieir origin can be obtained from the review by Grim. 



ITie basic building blocks of the clay niinermls, kaolinite, 
illites and mDntmorillonite consist of alteinate layers of gibbsite 
and eilica. Kaolijiite is a typical two layered clay containing one 
gibteite layer and one hydrated silica layer. Montmorillonite is 
conposed of a gibbsite layer between two hydrated silica layers i 
and is thus desigrsted as a three layer tvpe clay. 

2 . 2 Origin of Ion Exdhange Properties 

Soil particles are anphoteric as evidenced by their ability 
to bind both cations and anions, but, generally, they carry a net 
negative charge \^ich can be neasured by the zeta-potential technique. 
Ihe diarge is produced in two ways ; by isoinorpho^fi ion substitutions , 
and by ionization of hydroj^l groups attached to silicon of broken 
tetrahecfron planes; that is, 

HO + Si-OH^^ SiO" + FLO 

Ihe amphoteric nature of clay requires the existence of 
positive charges yiich may originate from hydrous oxides of iron, 
aluminum, and manganese and frcm exposed octahedral groups which 
react as bases by accepting protons from the surrounding soil solution, 
•ftus aqulring a positive electric charge. 

Because of the irregular shapes of the clay particles, the 
ciiarge density is higher on edges, comers, furrows and cavities than 
on plane surfaces. The charge density also varies with pH and the 
degree of ionization. 

2.3 Exchangeable Ions 

The electric diarge on soil particles is neutralized by 
exchangeable counter ions of an opposite charge and these are held 



i 



to the s:urface by 'Coulomb foarcee and Van der Waals -London forces. 
iQafmion exdiaaigeable cations are Ca,, Mg, H, Al, ' ' K, Na', and M . 
Calcium ions are generally predoiniiiant althou^ in, very acid soils 
A1,C0H)2_ constitutes a considerable proportion..^, incre-asing mth, 
deereas.lng pH., Sodium is usually veiy hi^ in alkaline soils. 

Ex.diangeable anions .are SO", Cl7 N0~, H^PO",, HCO" and , 
anions of humic acids. Some O'f tiiese ions dO' not always fijn,ction 
as counter' ions but are :mei:ely present in the soiil solution. In 
acid saline or sulphate soils the conte.nt of SOj^ is usually very 

Most studies on soil .anion exdnange have used PO^ and thus 
the .infoniation on -dhe exchange of this anion is more extensive. 
These investigations have indicated that the anion adsorptive powers 
of clay soils .are .assO'Ciated with the presence of Fe .and .Al conpounds 
sin,ce removal oif these .S'i±>s.tanoes by various lethods tends to either 
coiip.letely destroy or markedly red.uce anion adsorption ■( 9 ) . Th,e 
ledianism for anion .ads.orption. by P'Ure clay minerals proposes the 
'replacement O'f OH from, clay mineral surfaces,, especially around 
ttm flake edges, and adsorptio'n by anion-exdiiange sites existing 
.along the baS'Sl plane surfaces. ReplacQient o.f .adsorbed anions which 
have not been adsorbed, vi.a chemisorption is brought about by 'varying 
the pH. 'Ihis, reaction favors the theory of deactivation o.f tiie Fe and 
Al conpoimds present. 

Studies on phosphate fixation by soils and clay minerals 
indicate 'that "flie retention, of P0'~ is a result of the, reaction with 
oxides of Fe and .Al, the formatio'n of inso.l'uble salts of Fe, Al,, Ca 
as well .as fixation by the .clay minerals. The fi?atio'n by the clay 
mi,nerals is lo/j^ ranging from 0.03 to 0.07 millimoles P0~ per gram 
of clay,. Fixed P0~ ions can, be replaced or released from soil c,la,ys 
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by divalent arsenate, silicate, citrate, oxalate, hydrosiyl and 
fluoride ions. Citxmte and oxalate operate by dissolving iron 

eompoiinds , hydroxy 1 by decoinposijig clays , and arsenate and silicate 

( 8 ) 
by typical exdhange react iors. 

Because no single mechanism is responsible for "tiie uptake 

of anions in soils and clays, the temis ion exchange, adsorption 

and chemisorption , which are all applicable, should be replaced by 

the term "sorption". 

2.M- Ion Selectivities of Clays 

The selectivity of a clay refers to the preference for one 

counter ion over another. The valences of the counter ion have a 

. . . (4) 
strong effect on the selectivity and thus on equilibria. This 

effect is purely electrostatic and arises even in ideal systems. 

(4) 
As a rule, the ion exdianger prefers the counter ion of hi^er 

valence. The preference increases with dilution of the solution 

and is strongest with ion exchangers of high interfial molalil^. 

For anions, selectivity is in the order: 



POr > SO" > CI 
4 4 



and, for cations i 



Th > Nd > Ca > Na 



Exceptions to this rule are the anions in order of preference, 

SO^ > l" > NO3 > CrO^ > Br". 

For the case of ions of equal charge , the ion mth the siiB.ller 
hydrated radius is preferred. In tiie case of the alkali group 
the order of preference is, 



Cs* > Rb* > r > Na"*"' > Lis, 

■although the crystallo,graphic radius increases da«i the group; 
eimilarlyg the alkaline earths display the same trend; 

Ba > Sr > Ca > Mg: > Be: 



The order for monovalent anions is 



CMS' ' > ClOj^ >I" > NO'g > Br > 
HSO^ > NO'2 > Cl" > HCO° > 

OigCOO'" > 00' > F". 



3 . Sorption, Isotherms 

Sorption isotherms, describe, in the fonm of functions or 

dia,graiiis for a given (GOns:tanfl temperature, the dependence of the 
equilibrium concentration of the solute on the sorbent on the 
concentration in the external solution. It is preferable to relate 
the conc'entratio'ns in both phases - ion exchanger ,and solution ~ 
either tO' 'unit wei^t of the solvent in the phase or to unit volunie 
of the phase. 

Using this convention, the so^rption isotherm is either of 
two forrns: 

i = f(m^') ..... (1) 

or c = gCc. ) .......... (2) 

liiei^ m,. = concentration of solute i in millimoles per 
igpani solvent (molality); c = Goncentration of solute i in moles 
per litre (molarity); quantities with bars refer to sorbent phase. 

A single po'int on the sorption isotherm giving the ratio 
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of the ordinate value to 1±ie abscissa value is kno&m as the 
distribution Goefficient, It is defined as; 

<t ■' ■ '«i « • m » ii «' \ 'O / 



4 w *«•#»•«• \ H / 



yiere X. is the molal distmbution coefficient; 

1 

X. is the molar distribution coefficient. 

1 



3 . 1 Sorption Isotherro Models 

Early studies on adsorption equilibria were enpirical in 
nature, and attempted to fit eirpirical or s^ni-enpirical data to 
experimental results. The first major atteipt to give a quantitative 
forfrajlation to ion exdiange equilibria was Hade by Gans using 
the hhss action law in its sinplest form (ie: without activity 
coefficients). His relatiOTi was 



m 



'a"^ 



m 






const . 



»■'•«• « V' «i • m \t \J J 



■fliis equation holds for zeolites with regular bdiavior. Wiegner 

(12 13 14) 
and Jenny ' ' applied a form of the Freundlicti Isotherm, 

which contained two empirical constants, to ion exchange, nanely; 



m ~ I c^ - c 



1/p 



.*........( 6 ) 



where — = counter ions sorbed per unit wei^t of exdhanger; 



m 



c - equilibrium concentration of exchanging ion in solution | 
c - initial concentrHtion of exd^anging ion in solution | and 
k and p are empirical constants. 



1.1 



. , (15) 

A more convenient equation suggested by Vageler 

resenibles the Langmuir adsorption isotherm 



2L = Q Q 

m ^ Q + c' » (7) 



where Q = amount of exchangeable ion in solution; Q = weight capacity 
of exchanger; c' is an enperical constant. 

This equation is used in modified form as 



X _ ^ be ,,»,,,,.., ( 8 ) 
m ~ - 1 + be 



v^ere Q is the wei^t capacity of the exchanger; c is the equilibrium 
concentration and b is an empirical constant related to the energy of 
adsorption . 

For multiple ionic species , Boyd , Shubert and Adamson 
set up the follcwing equation by fornal analogy with the Langmuir 
theory: for species A of an electrolytic solution, 

/x\ ^Va ....(9) 



where — , = amoimt of A ion adsorbed per unit weight of exchanger; 

+ + , . * 

c., Cq = activities of cations in equilibriim solution; b. and b^ = 
constants related to the energy of adsorption of the cations; Q is 
the total exchange capacity. 
For the case involving ions of unequal charge, 

M "^aK ] ^ ..(10) 

m I A / v+iv + i "^ +VV+ 

where v and v are the valence states of the ions involved. 
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3 , 2 Phosphate Soiption Isotherms 

Most studies doiie on phosphate sarption isotherms have 
been done from the standpoint of ap*icultiiral evaluation. More 
specifically, their utilization has been directed tcward predicting 
the quantity of fertilizer phosphate required h^ a soil so that 

there is sufficient phosphate in the soil to produce a concentration 

(17 18 19) 
gradient necessary for net movement to the roots. ' ' 

The applicability of the Langmuir equation has been 

Cl8 19 20) 
strengthened by several studies , ' ' and a high positive 

correlation between the maximum value of P adsorbed and the clay 

content of the soil has been observed. 

The approach of the Langmuir analysis to soil adsorption 

readily allows one to evaluate the effectiveness of soils in limitini 

discharges from septic tank tile fields, Fur^tiier applications of 

eoil disposal are in the fields of slow sand filtration and land 

(22) 

disposal by spray xmgation. 

H , Epiamies of Contaminant Movement 

To predict the movenent of contaminants in soils the 
effects of dispersion and solution velocity must be evaluated to 
include the boundary conditions of the systeE, Many matheiatical 

node Is have been proposed for analyzing the dynamics of soil systems 

(21) 
with various reacting chemical species. The basis for these 

matheiTBtlcal descriptions is the equation of continuity for chenical 

species which is an unsteady state description of the system mass 

balance. 'Hie general form of ttie eqi^tion of change for species *'i" 



.^(2it) 

IS 



3c. 2 
- — + v*V c. - D. V c. + 1. (11) 



3t ^ -„ 1 1- 1 - 1 
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where c. is the concentration of species i; v is the superficial 

flow velocity vector; D. is the dispersion coefficient (replacing 

. . . . (25) \ 

the diffusion coefficient to include the mixing phenomena ) ; and 

R, is a rate of addition or removal of species. 

For dispersion of a single species of ion with no removal |. 

equation (11) heoomesi 



3c _ J. 3^c 



at 



9x' 



3— '- V 



9c 



8x 



(12) 



viaere c is tJie concentration of IJie ion; t is time; D is the longitudinal 
dispersion coefficient; v is the seepage velocity. For the semi-infinite 
boundary conditions, 



c(x5 0') =^0, X > Q 

c(0, t) = c , t > 

■' o — 

c(««5 t) = Og t > 



the analytical solution of the dispersion model is 



Q, 



1 - C 



2(Cti) 



^ "" '^"^(-ir) 



erfc 



1 + g 

2(Cn)' 



. (13) 

where the two diirensionleBS groups are 

1, -Oie Peclet nuntoer Pe = l/n = vL/D 

2. the pore volume C = vt/L 
and L is the leng^ of the model. 

If the solute is adsorbed on the medium, then the material balance 



3c 



dt 



= D 



d'c 



3x 



7 V 



9c 



3F 



9x 



(p 3t 

.:. ...(14.) 
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wher« i is the porosity of the medium; F is a form of the adsorption 

isotherm. 

The three possible adsorption isothems are as follows: 

1, Linear equilibriun ^ere F = k_c for whidh equation 



(14) changes to 



1 + 



/ ^l\ dc 
\ $ / 9t 



D 






9x 



m m '.H' "A m » • « m. \ JL<J' # 



» . (27) 

Ihe solution of Lapidus and Ammdson is 



- — = h erf c 

G 

o 



1 - C ' 



+ % exp { — — } erfc 
ri 



1 + V 



2(Cti) 

*«•«*•■««• V JL.U / 

where the modif ied pore volune , C , is e^ressed by 



V = 



(1 + k,/^) 



■m #'. 'i*j *.«*»«.« t \ . J^ / / 



fl2 



2, For Freuidlich equilibrium adsorption, F = k o ^^ 
a nonlinear parabolic differential equation resiilts| 



1 + J^^ c"2 - 



9 c 8 c 



Be 
3t 3x " ^ 3x 



■ « ■ . • ■ ■ • • • (18 ) 

A numerical procedure is required for the solution of this equation, 
3. Langmuir equilibriiM adsorption is given by 



(28) 



ae. 



P _ Qbc _ 

1 + be ' 1 + be 



ip ■ # tfe B « * 9 * « \ jL^' / 



ifiieTe a and b are constants. This ^xDduces* 



1 + 



Hi + bc)^_ 



3c 

at 



= D 



3c 
3t 



- V 



le 
9x 



^ ...(20) 
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Again 5 a namerical solution is required. 

% . 1 Bilinear Adsorption 

In addition to adsorption, theme is a continual exchange 
of ions. To include this ease properly ^ a kinetic model is required 

to account for tfhe bilinear rate of adsorption. Accoixling to Gupta 

(29) 
and Greenkom, this type of model is physically the meet appro- 
priate for adsorption in eoils and columns of soil. The adsorption 
reaction for chemical species, A, is represented by: 

A + eorbent -*- A'sorbent 
The rate of adsorption at soil stirface is 

(A)«(sorbent) j^- — ( A-sorbent) 

ad 



9t Kin 



] 



■..•«.••.,( 21) 

The concentratiQns are 

(A) = e 
(A*sorbent) = q 
(sorbent) - Q - q ,..,,..,..( 22 ) 
Substituting the rate expression in the raaterial balance yields 



(23) 



^ = K^J cCQ - q) - -4— q } 



3t jkin ^ H jr 

Ad 



with semi-infinite boundary conditions, 
c(x, 0) = 0; X >0 

q(x, 0) = Oj X >0 



( 7H 1 

'.• • « ■ • w 41 m m m- \ £.r^ J 
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c(0, t) = ci t >.0 
o — 

cC«, t) = 0; t >^0 
and K, , is the constant b, in tiie Langrauir equation. 

To obtain a solution, these equations must be solved nunErically on 

(29) 
a conputer. The solution of these equations is shown in Figijane 

1 as plots of c/c vs I/hs the nmriber of -pore volunes. These curves 
o 

shew the effect on breakthrou^ of tiie kinetic rate constant ^ the 
fluid velocity and the dispersion coefficient. 

Other inai^ematical approaches and conditions involving the 

tiransport of chemicals in porous media have been published by Lindstrann 

. ^vOO, 31) 
et al. 

Due to the difficulties in solving these equations, and 

the fact that a concentrm'tion profile is not necessarily required 

for field conditions, a simpler matiiCTnatioal technique is required 

to deterTOine breakttaxju^. 

H,2 Macroscopic M^s Balance 

In order to do a total mass balance for determining break- 
throu^ times for a soil-contaminant system, or, in o1±ier words, to 
deteimine the lifetime of a soil adsorption tile field, certain 
generalizing assumptions must be made based of tiie behavior of a 
conventional soil column. Cbnsidering Figure 2, a soil coluim is 
shown with a plot of the concentration-distance profile. The con- 
centration involved is that of the non-dimensional form c/c where 

o 

c is the average input concentration and c is the average concentration 
at a point x from the top of the coluim. For clarity, the averBging 
bars will be emitted in further analyses. Thus the following canr' 
ditions hold" 
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Figure 1, Effect on breakthrough curves of (a) kinetic 
rate constant (V = 0.05 cm/min, D = 0.05 cm^/min, Cq = 
0.00002 g mole/cm3, K^d = 540,000.0 cm^/g, and Q = 
0,00012 g mole/g adsorbent; for curves 1,2, and 3, %iji 
Is 0.002, and 2.000 cm3/g mole /min , respectively), 
(b) fluid velocity (D = 0.1 cm^/min, Kijin = 2.0 cm3/g 
mole/mln, Cq = 0.0002 g mole/cm3, lad = 540,000.0 cin3/g 
mole, and Q ^ 0.00012 g mole/g adsorbent), and (c) 
dispersion coefficient (V ^ 0.20 cm/min , K^i^ - 2.0 cm^/ 
g mole /min, Co ^ 0.0002 g mole cm3 , K^d ^ 540,000.0 qm? f 
g mole, and Q = 0.00012 g mole adsorbent). 



« « » i 



C/Co = 




€/Co= I 



G/C&^l 



^c/c« 



la) 





Flgyre 2. Adscrpliori Profiles in Soil Colymn 



•;] 



t - 0, X = 0, c/c = 1.0 
t = 0, X = h^, c/c = 0.0 

D O 

vAiere h„ is the bed depth., and 

B 

t >0, X = 0, c/c.^ =1.0 

t >0, X = x', c/c = (1...0, O.G) 

o 

ie." at some value, x* , c/c lies in the Pange fnom to 1.0 as 

described by the solution of the differential equations- Ihis range 

v^er^ c/c varies from a ■maximum, to its minimum over a bed incrcement 
o 

of Ax is kncwn as the reaction zone. 'This is the region where the 
actual adsorption and ion exchange processes are taking place. As 

the soil particles reach equilibrium with the inco'ming feed, this 
section of the bed becomes exhausted .and the reaction zone moves 
down the bed as depicted in Fi,gure' 2(b). After a period of time, 
the reaction profile will reach the bottom of the bed and the con- 
centration of the contaminant will rise from zero^ to a maxajnum as 
the concentration profile eierges from the bed indicating a depletion 
of the active sites in the soil (Figure 2(c)). 

To avoid differential equations and reduce the complexity 
of the bre.akthrO'Ugh phenoraena, the plug flow assumption must be used.. 
In this case the concentration .gradient is a step-function with a 
non-finite slope,, and thus,, dc/dx is iion-e,xi,stent since 3c/ Bx ■> «. 
The soil column behavior, thus,, is depicted in Figure: 3. 

4 . ,3 Derivation of Plug Flew Mass Balance 

Consider Figure 3(c) ijiiTiediately before bre.akthrough such 

that the d,imensionless concentration at the top of the columi is 



x. = c /c = 1 (25) 
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and at the bottom, the effluent is 

Xg = c^/c^ - ....(26) 

i.e. the column is e^diausted but bre.akthrough has not occurred. 
ITierefore, the col'uiinn mass input is V. c ovcv d .rixt.xi 
time period. At. 
The output is 

V ^c^ = V ^ X = ..., ...(27) 

out f out 

where V. and V ^ are. the total volumes in and out of the column 
in out 

over the tiire At, c is the input concentration and c^. is the efflue,nt 
concentration. For plug flow just before bre.akthrough , c^ = 0'. 

Imrae'diately on breakthrough , c_p = c .. 

The accumulation in the coluinn is 

where e is the porosity; ¥„ is the volume of the bed, ip„ is the density 

of the bed, a is the sorption equilibrium value of the contaminant on 
the soil. 'The first term of the acciinulation is th,c amount of cxm- 
taminant contained 'in the voids of the bed, and tlie stxx^nd ttmii ir> Iht"* 
anrount .adsorted on the soil at breakthrough. 
Equating, 

^in^o= "Vo" ^'Pb^^B (28) 

V . c - eV„;C = q ip.eVn ... . ( 29 ) 

m, o Bo ^«> B B 



m^^ 



V^-eV =W^ C30) 
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O.' 



Dividing by eV„; 
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If the cross-sectional area for flow is A over a tljne 
period At where A is the cross-sectional area of the coluiiiin, the 
flow velocity of the aqueous phase is 



V = V. /e Mt (33) 

m m 



The rate of movement of the breakthrough front is 

V„ h„ 
V _ B _ B 

^ Mt ' M ..(34) 



Therefore dividing numerator and denominator of equation (32) 
by eAAt, ie. maltiplying by unity. 



1/eMt _ ^ 
1/eMt 



yields 



"= JL ^ ■ " ■ ' «, "i; ^.' <i >i( i» -.*■; -w m • m ^,00/ 



V EC 

C o 

which relates the flow velocity throu^ the bed to the velocity 

of the moving front. 

The quantity, q^p„/c , is a dijiensionless coefficient which defines 

the liiniting distribution between the solution and sorbent phases. 

If q is replaced by the adsorption isotherm 3Ji the form 

Obc 
X _ ^ o 

^ " m " 1 + be ......... .(8) 

o 



?b 



equation (35) becoies 

V Qbc p 



u 1, ^ ^ o ■ B ,....(36) 

— ^ _ JL + — — 

C 1 + be EC 



\l = 1 + Q^'p 



B^ / 1 1 ......... ..(37) 



^c E \1 + be 

Two oas;es may result depending on the nia,gnitude of be : 

(i) be « 1 

o^ 

\ ., , Q^Pb (38) 

V 

Q e 

and, thus the velocity ratiO' becomes independent- of the influent 

conoentration , 

(ii) be » 1 

o 

^ = 1 + B .... (39) 

c ec 

o 

for which the velocity ratio depends on the 'neximuB ads.orpti've' 
capacity Q .and the influent concentration,,, .and not on the isothenn,, 
■parameter "b". 
5 . Qi Application of ContaiTiin.ant Movgnent Equation, to Field Conditions 

Ej^erimental data on flow of a contaminant from a source 
have indicated a dispersion fo^r unconfined flow which is basically 
linear ,in na,tu]ne. The extent of dispers.ion or the dis,persion co- 
efficient has been shown to be a function of the flow Reynolds number., 

(23) 
Ihe dispersion and the Reynolds nunijier have been correlated xn, the form, 

D = a^vRe^^l ...... (i+O) 

where D is the dispe^rsion coefficient, .a and b^ are constants, v is 

the kinematic viscO'Sity and Re is the Reynolds nijmber based on pei''-'- 

meability. In ino:s,t cases b^ h.a:r. tlic value of .approximately unity. 

Consider Figure 4 in which a source of pollution is 

entering the soil at a concentration c through, a cylindrical area 
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Figure 4. Contomirartt Dispersion Model 
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at point 1 and then disperses so that at seme point at a distance, 

r^, fron the source, the flew area has expanded as shown by point 2 

These areas subtend angles and a whose vertices are located at a 

distances , r and r ' , back from the source . 
o o 

At point 1 the arc length and depth are 

a^ = r 9 : K = r' tan a , (41) 

1 o 1 o 

and at point 2 

^2 " '"o ■ "'" ' ^'2 '"o 



a^ = (r_ + r)0 ; h^ = (t[ + r)6 tan a 



(H2) 

Therefore, the flow areas are: 

point (1) I At = a-,h-, = r r' tan a 

X L L O O 

*,........ (43) 

point (2): A^ = a„h„ = (r + r) (r' + r)0 tan a 
- 2 2 2 o o 

(44) 

The relationship between concentration and flew areas is 
given by -tiie integrated continuity equation: 

^A^l" ^2^2^2 -^^^ 

for two positions, 1 and 2. 

, • . c^ = ^x^Vj^/A2^2 ' (46 ) 

For simplification purposes , it must be assuined that 

v^ Z v^ (47) 

neglecting differences in velocity resulting from lateral dispersion*, 
there forB , 

Substituting values for A., and A from equations (43) and (44) 
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or 



c^r r'0tan a 
loo 

(r + r) (r* + r)0tan a 
o o 



c_r r' 
loo 

(r + r) (r' + r) 
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which gives the concentrHtion at a distance r from the source. 
Applying this result to equation (39), the following is obtained; 



V 



^ = I + £_ 

V e 

G 



1 + be 



o 



oir + r/l r + rj 
\ o / \ o / 



This equation includes the effect of distance on the rate of 

contammant inovement. The values of r and r' are given by 

o o *^ -^ 



r = 
o 



W/2 



r' = 
o 



sin©/ 2 "' tan s 

*.. (52) 

and h^ is usually equal to the tile diameter as a first 

approxifiHtion, a is an experimentally detemiined angle, W is the 

width of the tile bed and is the plan view angle dependent on 

dispersion. 

If dispersion depth values h, and h^ are leicwn. 
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(53) 



For a constant depth above an impenneable layer, h, = h„, 1±ien 
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In m\y case, some approxiination to the contaminant flow profile 
is required to be made based on the flow hydrology. 

6.0 Experimental 

6.1 Detemnination of Phosphorus Existing in Soil 

'The concentration of phosphorus existing on the soil of 
a tile field varies with distance from the point where it initially 
leaves the tile to enter the soil. To observe the gradation in 
phosphorus concentration with distance, soil samples taken at 
several points away from the source of contamination miist be analyzed 
to determine the amount of inorganic P present on the soil. 

Most soil phosphate detemiinations involve two distinct 
steps; extraction of the soil to obtain the phosphorus in solution, 
and tiien a quantitative analysis of the phosphorus content. Colori-^ 
metric raetliods forming the mol^todenun blue complex appear to be the 
most sensitive technique for soil extracts containing snail amounts 
of POj^ - P. 

6.1.1 Principles of the Extraction Technique 

Soil extraction using a strong acid such as IM sulphui''ic 

acid gives an approximation of the total inorganic phosphorus in 

(32) (33) 
a soil. Parton suggested that the solubility of inorganic 

phosphorus is a function of pH and the total amount can be brought 

into solution provided the pH is made sufficiently low. At pH 1.0 

all the inorganic phosphorus will be mobilized but not necessarily 

(33) 
in solution as soiib readsorption of phosphorus will occur. 

Parton overcame this problem using a radioactive tracer to determine 

the degree of readsorption. Because the time to reach equilibrium 
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may be in the order of three weeks, the effects of r>eadsorption 
have been neglected in this study. 

Ihe phosphates in solution are generally considered to 
be in the fomi of aluirdnijii phosphate, iron phosphate, and calcium 
and dicalciiin phosphate. These phosphates may be^ frmctionated 

Beparately with different extractante although there is a gCTieral 

(32) 
degree of overlap of each phosphate form. 



6,1,1 Reagents 

IM H,,SO,, (56.6 ml 18M H„SO,, (cone) in 1 litre). Reagents 
according to "Standard Methods" for POr' - P analysis. 



6*1.3 Procedure 

Wei^ out accio^ately 1.0 pn of soil passing a 60 inesh sieve 
into a 250 ml erlerm^er flask. Add SO ml of IM H„SO^ solution, seal 
the flask and shake for h hour. Because of the presence of sewage 
on the soil, H^S may form, and must be released before pressure builds 
up. Filter the solution conpletely throu^ a O.i+5 y millipore filter. 
Mash the soil residue twice with 10 ml of distilled water. Dilute 
the entire sample to exactly 100 ml in a volumetric flask. Analyze 

itti aliquot for inorganic P usijig the stannous chloride tecJmique 

.■ . ,- (31) 

described in "Standard Methods", 



6.2 Phosphate Isotlierm Determination 

The previously described phosphate analysis measures the 
amount of phosphorus or phosphate already held in the soil. The 
isotherm test measures the aniDunt of phosphorus a soil will rapidly 
take up# 



il. 



It is also worthwhile to note that the maximum capacity 
for phosphate adsorption is not equivalent to the soil anion ex- 
change capacity. It is possible that the sim ot the phosphate 
ions present on the soil .and the adso'rption capacity ni-iy approac^-^h 
the value of the anion exchange capacity provided that the amDunt 
of adsorted chloride,, sulphate and nitrate is low, 

6.2.1 Reagents 

Standardized P solutions of concentrations 1.63 mg P/1 

C5 ing/l POr) made with deionized water, and deionized water containini 

OuO'lM NaCl or O.OIM CaCl . Reagents for the stannous chloride 

z 

(34) 

technique are listed in "Standard Methods". 

6.2.2 Procedure 

Wei^ out accurately, approximately 1, 3, 5, 7 and IQi pn 
sanples of soil passing 200 mesh (or 60 mesh in the case of sand), 
into 250 inl erlenirp/er flasks. Add 150i ml of standardized phosphorus 
reagent, seal the flasks, and shake for 24 hours (as a standard 
shaking tiiiB) or longer if desired. Allow to stand 1 hour, then 
centrifuge the supernatant to obtain lOiQ ml O'f solution. Filter 
the centrate and analyze an aliquot for ¥0. - P using the stannous 
^chloride method. 

As mentioned in section 6.2.1 standards containing 0.01 M 
Ma CI or CaCl« were used because they siimlate soil solutions, 
stabilize pH, .and aid in the coagulation oif soil colloi.ds which could 

Interfer^^ with phosphate nneasurements . '" '■ However it lias been 

(36 37) 
observed by several authors ' ' " that calcium contacting 

solutions gave hi^er adso-rption values than sodi'um chloride 

contacting solutions, especially at pH values ,^eater than 5.5. 

'BiUB,, the NaCl solutions are usually favoured. 
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7.0 Results and Discussion 

7.1 Phosphorus Analysis' of Soils 

To test the analytical technique of determining soil 
phosphoruS' , soil saiTples were removed from nine points in a 43.2 
cm colunm of sand throu^ which sewage had been, flowing for 
one year. The phosphorus in the feed varied .from 1.5 to 35 mg/l. 

ITie sanples were analyzed for total inorganic phosphorus 
as de.s.cribed in section 1.4,, and the results plotted vs. distance 
as in Figure 5. The curve drawn throug|i the points shows the 
approximate phosphorus profile in the column. The shape of the 

sorted phase concentration curve is similar to that predicted by 

(30) - . . . 

Lindstrom and Stone ' jji their mathematiical .analysis of this 

type of probleii. Their curves for various exchan,.ge •coefficients, 
•are rep'roduced in Figure 6^ for coiiTpar.ison . 

'The existing background concentration on the soil may 
be considered tO' be 105 pg,/giTi whi'Ch is the bottai point of Figiffe 5. 

In a similar' manmer, .s^aonnples from a septic tank tile field 
.may be analyzed to detennine P0^ - P profiles and thus deterTiiine 
the 'e-xtent of utilization of the soil by the septic tank effluent. 

A study of this nature carried, out over a long period o.n 
a specific tile field will also give information on the rate of 
mDvonent of the contaminant front, and thus predicticns as to the 
life of the tile field may be made. 

7.2 I%osphonis Isotherms 

7 . 2 , 1 Langmiir Isotherm Model 

A plot of phosphorus adsorbed and equilibrliin concentration 

of phosphorus in the supernatant, in units of i^g, P/lOO pn soil and 
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Figure 5. Phosphorus Concentrations in a 43.2 cm 

Laboratory Column 
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Figure 6. Relative sorbed phase concentration 
S ( X , 1 )/')rUQ distributions for some values of K 
(rate of exchange coefficient) as a function of 
distance x (cm) from the top boundary (x = plane) 
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mg P/1, respectively, is shown in Figure 7. It may be observed 
that this plot of x/m vs c produces a curve that approaches a 
horizontal as:3mptote. Considering the Langnuir equation, 



Qbc 



m 1 + be 
e 



it is obvious that as c_ becomes large, 

ie be »1, 

e 



^ .,,.. (8) 



m 

indicating that the isotherm curve will approach the value of Q, 
vdiich is the maximum sorptive capacity. 

To determine Q more: 'readily, the Lan,gmijir equation may 
be expressed as 



1 + be 
e 



-^ Qbc 

111 ^ e 



.(55) 



— (56) 



Qb \ ^e / Q 

11 ■. . 

Plotting vs — — generally produces a straight line the 

IT 

intercept of which on the l/(x/m) axis gives 1/Q, allowing one to 

obtain a value for the maxirr^'um adsorption capacity , Q , ( Fi,gure 8 ) . 

The Q value is frequently cited in the soils literature as a 
^ max ^ J 

characteristic of the soil; however, the x/m value of most practical 
interest is the one that occurs at the concentration of the waste- 
water to be applied in the field. 

Other examples of isotherm plots are shown in Figure 9 for 
various soils from the Kingston Ontario area. The variation in Q 
max is due to the different compositions of the adsorbing fractions, 
which consist mainly of clay and silt, and the amount of phosphorus 
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Rgure 7 . Phosphorus Adsorption Isotherm for Soil 
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O 73% passing 200 mesh 
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Figure 9. Riots of i /(xitn) vs l/Ce for several Kingston, Ontario soils 
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previously adsorbed and occupies active adsorption or exchange sites. 

7.2.2 Active Adsorbing Fraction of Soil 

(38) 
According to Bear it is well kncwn from numerous 

investigations that the ion exchange property of soil is due almost 

entirely to the clay and silt fractions (<20 ijm) and the organic 

matter; the colloidal material of soil beir^ most iinportant. 

Figure 10 shows a plot of Q max vs the percent silt for 

various samples taken from two different areas. In both cases 

there is a correlation between % silt and clay and the value of 

Q^ _ for phosphorus . Studies done on sand of high silica content 

showed no phosphorus retention. Also soil fractions larger than 

100 mesh containing no organics exhibited negligible phosphorus 

removal. 

7.2.3 Phosphorus Desorption 

If physical adsorption were the dominant mechanism, adsorbed 
phosphorus would be rapidly desorbed when the concentration in the 
applied solution is low. Samples which had undergone the adsorption 
process were subject to desorption using distilled water and shaking 
over period of one week. These samples did not release phosphorus 
of an equivalent amount as adsorbed but retmned a significant 
quantity. This is the criterion which indicates that the removal 
phoiomenon is not solely the result of physical adsorption but 
other mechanisms also remove or insolijfclize the phosphate. 

Rainfall leaching will tend to remove a fraction of the 
sorbed P, but this effect is usually minimal since rainfall per- 
colation rates are usually quite small relative to the wastewater 
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application rate, and, witfi the long timie involved in .ground water 

travel, considerable ■retention and slow mineralization of phos:phate 
ion occup, 

7.2,'^ pH Effects 

The overall effects of soil pH on adsorption are^ well 

icnoWTi and have been thorougily discussed. '"'' In general, the 

reioval of phosphate decreases with, inCT^asing pH except for the 

conditions involving a CaCl^ solution at hi^ pH. The literature 

confimis this trend for acid soils in which iron .and aluJidnum are 

responsible for tiie removal. ' "' ' '"' - However for alkaline soils 

in whi'Ch calcium is responsible for ■renovals , the reverse is true 
of ph/ 39, HO) 

The data from the Kingston soil samples ,are plotted in 
Figure 11 shewing the trend of x/m with soil supematant pH, As, 
anticipated with these acid soils, the .amoimt adsorbed decreased 
with increasing pH. In general, no prediction of this variation 
could be :made since the independent variables such as initial 
soil pH, initial phosphorus concentration, and the nature of the 
soil counter ions and co-ions are dependent on the nature of the 
clay and silt mineral content. 

7.2. 5 Effects of Contaoting Solutions 
Isotheim analyses were performed on a To'^ronto soil of 
pH 6.8 ,and passing 200 mesh, using different cantacting SDlutions 

of phosphorus. Ihe solutions used were: 

1. phos^orus in distilled water 

2. phosphO'rus in OuOlM NaCl 

3. phosphO'rus in O'.OIM CaCl^ 
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Figure It. Variation of (x/m) with pH 
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As mentioned before ,, sciie investigators favour using 0,0 IM NaCl 
or O.OM CaCl„ because it nmy simulate soil solutions, stabilize 
pH,, and aid in coagulation of soil colloids which could interfere-: 
witili, phosphorus iieasursnents « 

Ihe isotherms for the -diree contacting solutions are 
plotted in Figure 12 . The valuee of 5^ for this soil are sho^^m 
in Table 1. 

Table 1. Effect of Contacting Solutions 

Solution imx 

1. P in distilled water 10.3 

2. P in O.OIM NaCl 26.1 

3. P in O.OIM CaClg 28.2 

Ihe results show that Q is increased by the use of 

NaCl and CaCl„. The question of which solution should be used must 
be based on the. nature of the effluent entering the soil and the 

n ++ + 

concCTitrmtlon of Ca and Na already present m the ground water. 

+ ++ 

In "tile large, Na is preferred over Ca because of the precipitating 

++ — 

effect of 'Ca' in p'rodueing hi^er () , values.. In this investigation, a 

only phosphorus and 'distilled water solutions were used in order to 

minijTiize effects which were, not native to the soil s.ainple being 

analyzed. 

7.2.6 Approach to Equilibriijn 

TO' test the apprDach tO' equilibriun,, twO' sets of samples 
were weighed out and shaken for five days with one sample from each 
set ranoved and analyzed after 4 hrs, 24 hrs, 72 hrs and 120 hJB . 
One set of sanples contained 1.0 gram quantities of Toronto soil. 
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O PO4 in wotor + 0.01m CaCl2 28.2 

# PO4 in water alone 10.3 

A PO4 In wafer +0.01 NaCi 26.1 

Sol) pH = 6.8 
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Figuri 12 . Effect of Stobiliiing Solytions on PO4- P Adsorption 
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pB- :§..&,. .and' the o^'ther set contained 5.0! gram eamples of the same 

soil. The phosphorus solutions; oontained 1.63 mg/l P (5 mg/1 P0~) 
for both sets. Hie values for x/m and tijne are plotted in Figure 13 
for each set of samples. 

From the curves it is apparent that equilibrium was reached 
in 30' hrs for the 5 pi saonple, and almost approached at 120 hours 
foc" "ttie 1 ^am s^airple . One reason for ttie differences in appiroaidi 
to equilibrium is ttie variation in relative initial ratios of phosphorus 
to soil weights as shcsAjn in Table 2. 

Table 2. Phosphorus to Soil Ratios 

Sample Size Solution Initial ratio 

Cgm) Concentration CiBg P/1) P to: Soil (mg F/lQiQ gm.) 

1.000 1.63 22.3 

5.00.0' 1.63 4.84 

The idif ference in the rate of approa'Ch to. equilibrium is 
goverried by Fidk's First Law which is expressed as 



--— = k (S* - S) .,.,». .,...(57) 
at 11 



where S is the amount of P adsorbed per unit wei^t of soil, S^'* 

is the. equilibriiin amo^unt adsorbed and k is the mass transfer 
^ m 

coefficient for a particular batchwise shaking process. If we 
consider the 1cm range of the equilibrium isothenii, which is 
,appro:'3djiiately linear, the relationship between the solution con- 
iOCTftration .and the adsorbed concentration is 

S* = \c ^ **.*».•,.... .(Si) 
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v^ere X is a partition coefficient and c is the solution concentration 
■pef erred to the unit weight of saiiple. Therefore Fick's law becoroes 

^ = k ( Xc - S ) (59) 

dt m 

The value of c is 22,3 ing P/100 gm for the 1 gm sanple 
and 4.84 mg P/100 gm for the 5 gm sanple. If the initial amount 
of P on the soil is S , the concentrntion driving forces for the 
two sanples are ( X x 22,3 - S ) and ( X x 4.84 - S ), respectively. 
The 1 gram sanple has a greater slope, dS/dt, and thus, the time 
to reach equilibrium will be greater. 

This result thai irnplies that the equilibrium tijas for 
each sanple wei^t or concentration will be different, and attainihent 
of equilibrixjm is only reached after 120 hours for all cases covered 
in this concentiHtion range. 

For the range of concentrations involved in this stud^y, 
78% approach to equilibrium can be attained within 24 hours. The 
rate of water travel throu^ soil is often less than several meters 
per day so adequate contact time per day is available for P removal. 

7.3 Example Calculation of Contaminant Movement in Soil 

The application of the contaminant movement equation derived 

in Section 5.0 is illustrated below using a hypothetical example. 

Data for an actual situation are in the collection phase and are not 

yet available. 

Consider the case shown in Figure 14 for a septic tank tile 

field wi-|± 4 laterals 15.24 m in length and 1.829 m apart. The 

following physical parameters would be obtained by measurement: 

V = water movement rate = 0.1 m/day 
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Figurt 14. Tile Field Configuration 
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p^ = soil density = 2.7 pn/cm or 2.7 x 10 ^/l 

e = soil porosity = Q,2 

r = movement distance to body of water 

= 7.32 + 15.24 = 22.56 m 

The chemical parameters are: 

c = the septic tank phosphorus effluent concentration 

obtained by chemical analysis, 

5 - 4,52 mg P/100 gm or 4,52 X lO"^ mg P/gm soil 
b = 2.193 1/mg 
and C) and b are obtained from a Langpair isotherm analysis. 

A further assumption to be made is that the phosphorus 
penetrates to a depth of 61.0 cm at the lateral farthest from the 
lake, and increases linearly to a depth of 183.0 cm at the bank. 



ie: h, = 0.61 m 
h^ = 1.83 m 
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and the angle of dispersion, 0, is 28 . 

Fran the appropriate equation for the relative rates of movement, 



the value of v may be obtained. 



and 



V 



0) 



= 1 + 



^^B^ 



1 



1 + be 



r = 
o 



r = 



W/2 



sin0/2 

15.24/2 
sin 14^ 



= 31.5 m 



h. 



ol r + r I ho 
I o / 2 



(53) 



V 



U) 



V 



= 1 + 



H . 52xl0~^x2 « 7x10^x2 . 193 

0.2 



1 + 2.193x20.0 / 31.5 \ 0.61 
( 51+. 1 ) 1.83 



1 + 1338.2 (0.1051) 
141.7 



* , for V =0,1 m/day , 



.-H 



V = 0.1/141.7 = 7.058 X 10 m/day 
To move 22.56 m, the time required is 

. 22.56 m o nnr- ^ Tn^ J 

t = j7 — = 3.196 X 10 days 

7.058x10 

or 87.6 years. 

If the rate of water movement were 0.3 m/day, the time 

for the soil to be exhausted would be 29,2 years. 
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Conclusions 

1. Riosphorus removal fr>cm effluents passing through soil is caused 
by a combination of three mechanisms: ion exchange, physical 
aefeorption and precipitation. 

2. The ability to remove anions from solution is related to the 
structure and mineral content of the clays, 

3. Clays act in the ion exchange mode as conpared to the adsorption mode, 
prefer the counter ion of higher valence, and selectivity increases 

.. with decreasing hydrated radius. 

4. Phosphate isotherms can be determined by a 24 hour batch equi- 
libration with a solution of a kncwn phosphorus concentration. 

5. TTie Langmuir isotherm provides a convenient method for modelling 
the sorption variation with concentration. 

§■:- The langnuir model can, be readily extended to multi-species 

solutions to estinate sorption properties. 
7. Ilie rate of movonent of a contaminant in groundwater flew can 

be estimated using a form of the equation 

V /v = 1 + "^^B 
uj c ec 
o 

and "die required parameters are the maximum adsorption capacity, 

tiie soil density and porosity, and the concentration of the 

contaminant . 
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Nomenclature 

a enpiricalcbnstant = Qb 

^1 enpirical constant 

... 2 
A, cross -sectional area at position, i, m 



3- 



€3npirical constant related to energy of adsorption, l/mg 



b«. value of b for species i 

h^ enpirical constant 

e concentration, mg/1 or moles /I 

■Oi concentration of species i, rag/1, moles/l 

e equilibrium concentration, rag/1, moles /I 

c • 

3 final concentration, mg/1, moles /I 

c initial concentration, mg/l, inoles/l 

c' empirical constant 

D dispersion coefficient, era/sec 

h_ height of bed, m 

k empirical constant 

k^ linear equilibrium constant, l/ng 

It^ kinetic rate constant 

k mass transfer coefficient, sec 
11 



equilibrium coefficient based on the mass action law, J 

■ 1 

i 

lii length of bed, m -' % 

m weight of soil sample, gm 

% molality of species, i, m mDl/gn solvent | 

& molality of species, i, on exdTanger, ra mol/gm solvent | 

Tis depth of dispersion at point i, m |l 

n^i constant for Freundlich Isotherm 

p enpirical constant 

Pi .Peclet Number = vL/d 
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q adsorbed concentration = x/m, iiig/,gm; or (x/m)=^Q|Q, mg/lOO ,gm 

Q wei^t capacity of exchanger- mg/pn or mg/lOO 311 

Ip-' radial distance, m 

Re Reynolds' Number = vdp/iJ 

■R. rate of addition or removal of species i, mg/l sec 

S> anount of adsorption on soil, :irg/,pi 

S* aiTO'imt of adso'rption on soil at equilibrium,, ]mg/,gii 

% time, sec 

w veloci'ty, m/sec 

3 

f voluire, m, 

W width of initial dispersion, m 

tt amount of solute adsorbed, mg 

Greek Letters 

<m .angle, rad 

i porosity 

1 angle, rad 

1 1/Pe 

n. ■molal distribution coefficient 
1 

nl ndlar' distribution coefficient 

p vieeosity, cp 

m va:ience state 

I pore ■volijiie = vt/L 

5' nodified pore volume 

p: density 

p soil density 
B 

■§ porosity 
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